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ABSTRACT 

We present Atacama Large Millimeter/submillimeter Array (ALMA) obser- 
vations of the [C II] 158 /xm fine structure line and dust continuum emission from 
the host galaxies of five redshift 6 quasars. We also report complementary obser- 
vations of 250 GHz dust continuum and CO (6-5) line emission from the z=6.00 
quasar SDSS J231038. 88+185519. 7. The ALMA observations were carried out 
in the extended array at O.T' resolution. We have detected the line and dust 
continuum in all five objects. The derived [C II] line luminosities are 1.6x10^ to 
8.8x10^ L0 and the [C II]-to-FIR luminosity ratios are 3.0 — 5.6 x 10~^, which is 
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comparable to the values found in other high-redshift quasar-starburst systems 
and local ultra-luminous infrared galaxies. The sources are marginally resolved 
and the intrinsic source sizes (major axis FWHM) are constrained to be 0.3^' to 
0.6'' (i.e., 1.7 to 3.5 kpc) for the [C II] line emission and 0.2'' to 0.4" (i.e., 1.2 
to 2.3 kpc) for the continuum. These measurements indicate that there is vig- 
orous star formation over the central few kpc in the quasar host galaxies. The 
ALMA observations also constrain the dynamical properties of the atomic gas in 
the starburst nuclei. The intensity-weighted velocity maps of three sources show 
clear velocity gradients. Such velocity gradients are consistent with a rotating, 
gravitationally bound gas component, although they are not uniquely interpreted 
as such. Under the simplifying assumption of rotation, the implied dynamical 
masses within the [C II]-emitting regions are of order 10^° to 10^^ Mq. Given 
these estimates, the mass ratios between the SMBHs and the spheroidal bulge 
are an order of magnitude higher than the mean value found in local spheroidal 
galaxies, which is in agreement with results from previous CO observations of 
high redshift quasars. 

Subject headings: galaxies: starburst — galaxies: evolution — galaxies: high- 
redshift — quasars: general — submillimeter: galaxies 



Introduction 



Quasars at redshift 6 and higher provide a unique sample to study the formation of 
the first supermassive black holes (SMBHs) and their host galaxies at the epoch of cosmic 
reionization. There are currently more than sixty quasars known at z~6, sel ected froni lar^e 
optical and near-infrar ed surveys (e.g. the Sloan Digital Sky Survey [SD SSj.lFan et al.[ | 2006 |: 



Jiang et al.l[2008l . l2009 l: the Canada- France High Redshi ft Quasar Surve y, 



Willott et al 



20091 I2OI0I . the UKIDSS Large Area Survey [ULAS]. iMOTtlock et all EoOqI. I2OII 



Pano ramic Survey Telescope and Rapid Response System [Pan-STARRS], 



20071 . 



and the 



Morganson et al. 



(Kurk et al. 


2007; 


Jiane et al. 


2007; 


Willott et al. 


2010; 



2OIII ). The formation of 10^ M© SMBHs at the highest redshift suggests fast bl ack hole 



accretion and significant SMBH-galaxy evolution w ithin 1 Gyr after the Big Bang ( ILi et al. 
2OO7I : iNaravanan et al.lboosi : IPi Matteo et al.ll2012h . 



Dust and high-order molecular CO line emission in the host galaxies of the z^6 quasars 
have been searched fo r using the Max Planck Millimeter Bolometer Array (MAMBO) on the 
IRAM-30m telescope feertoldi et aPboOSal : IPetric et al.lboosi : IWang et al.lboOTl . boosl . boilal : 
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Qmont et al.l l2013r> and the IRAM Plateau de Bure interfero meter (PdBI, iBertoldi et al. 
2003bl : IWalter et al.l[2003l : ICarilli et aPboOTl : Iwang et al.lboiol ). respectively. The MAMBO 
survey of SDSS z^6 quasars at 250 GHz (typical la errors of 0.6 mJy) found that about 30% 
of them show strong continuum emission from dust at a temperature of 40 to 60 K, with FIR 



luminosities of a few 10^^ to 10^^ and infer red dust masses of a few 10^ M^t) (IBertoldi et al. 







2003al : IPetric et al.l l2003l : I Wang et al.l l2011al ) . Most of the FIR- luminous z^6 quasars are 
also detected in CO line emission, implying that there is ^ 1 ^^ Mc?) of molecular gas i n 
the quasar host galaxies d Walter et alfbood : ICarilh et all hopj . IWang et al.lboiol . boilbh . 
The strong FIR continuum and molecular CO line emission suggest that there is active star 
formation in the host galaxies of these FIR luminous z^6 quasars. The star formation rates 
i nferred from the FIR luminosities ran^e from a bout a few hundred to one thousand Mq yr~^ 
feertoldi et al.ll2003al : IWang et al.lboosl . boilah . 



The CO detections from the FIR-luminous z^6 quasars also provide constraints on the 



(Walter et al. 


2004; 


Wane et al. 


2011b) 



luminous quasars, SDSS J114816. 64+5251 50. 3 at z=6.42 (hereafter J1148+5251), u sing the 



Very Large Array (VLA) and the PdBI f lWalter et al.l l2004l : iRiechers et al.l 120091 ) at sub 



arcsecond resolution reveals a CO source size of about 3.6 kpcxl.4 kpc (FWHM) in the 
central region of the quasar host galaxies. Recent JVLA observations of the CO (2-1) line 
emission from another FIR-luminous quasar, SDSS J092721. 82+ 200123. 7 at z=5.77 . at 2'' 
resolution has constrained the source size to be about 10 kpc ( Wang et al. 2011bl ). The 



line widths of the CO-detected z^6 quasars range widely from 160 to 860 kms~^. The host 
galaxy dynamical masses (within the CO emitting region) estimated from the CO line widths 
and the sizes of the bright detections indicate a medi an SMBH-bulRe mass ratio about on e 



order of magnitude higher than the present-day value ( IWalter et al.ll2004l : IWang et al.ll2010l ). 



The PdBI image of [C II] 158 /xm fine structure line em ission from J 1148+5251 also sug- 
gests intense star formation and quasar-galaxy co-evolution ( Maiolino et al. 112005 : Walter et al 
20091 : iMaiolino et al.ll2012l ). The [C II] line from this object shows two components. The 
narrow-line component with a line width of 345kms~^ and a spatial extensi on of FWHM 



1.5^^ (8 kpc) traces the distributed star formation in the quasar host galaxy ( IMaiolino et al. 
2OI2I ) and the compact core of this component within the central ^ 0.3^' reg ion suggest a 
very high peak star formation rate surface density of ^ 1000 M© yr~^ kpc~^ ( IWalter et al. 



2OO9I ) . The broad-line component shows a line width of FWHM^2000 km s ^, distributed 



over ^16 kpc, suggesting gas outfiow driven by the central luminous quasar ( IMaiolino et al. 



20121; IValiante et al.ll2012l ). The detection of dust continuum, CO, and [C II] line emission 
in the quasar host galaxies at z^6 suggest an early phase of SMBH-galaxy evolution, in 
which the SMBHs are accreting at their Eddington limit and the quasar stellar bulges are 
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still accumulating their mass via massive star formation. The interferometer observations of 
J 1148+5251 have demonstrated that the [C II] hne emission is an ideal tracer of star- forming 
activity and gas dynamics in the nuclear region of the starburst quasar host galaxies at the 
highest redshift. 

The Atacama Large Millimeter/submillimeter Array (ALMA) in the early science phase 
provides a wide frequency coverage for observations of fine structure line emission at high 
redshift. It also provides the required sensitivity and allows detections of line emission with 
peak fiux density of a few mJy in only one hour of observing time. Additionally, the spatial 
resolution of ^0.7^' in Band 6 and 7 in the extended configuration can measure or tightly 
constrain the distribution of the gas component from high redshift quasar host galaxies. 
Thus, the Cycle phase of ALMA opens a unique opportunity to expand the observation of [C 
II] line emission to a large sample of quasar-starburst systems at the highest redshift. In this 
paper, we present new ALMA observations of the [C II] line emission from five FIR luminous 
quasars at z^6. We also report IRAM observations of the 250 GHz dust continuum and CO 
(6-5) line emission from the quasar SDSS J231038.88+185519.7 at z=6.00 (Fan et al. 2013 
in prep.). We describe the observations in Section 2, present the results in Section 3, discuss 
the star forming activity and gas dynamics based on the [C II] line detections in Section 4, 
and summarize the results in Section 5. A A CDM cosmolog y with Hp = 71km s~^ Mpc~^, 



= 0.27 and Q\ = 0.73 is adopted throughout this paper ( ISpergel et al.l 120071 ). 



Observations 



2.1. IRAM observations of SDSS J231038.88+185519.7 



The quasar SDSS J231038. 88+185519. 7 (hereafter J2310+1855) is a broad absorption 
line quasar discovered in the SDSS (Fan et al. 2013, in prep.). The redshift measured with 
the quasar rest-frame UV line emission is z=6.00±0.03. It is one of the optically brightest 
sources among the known z^6 quasars with rest frame 1450 A magnitude of mi45o = 19.3. We 
observed the 250 GHz dust continuu m from this object with the MAMBO-II 117-element 
array on the IRAM 30-m telescope ( IKreysa et al.l 119981 ) in the winter of 2010-2011. We 
adopted the standard on-off photometry mode with a chopping rate of 2 Hz by 32^' in 
azimuth. The total on- source observing time was 100 minutes. The data was reduced with 
the MOPSIC pipeline ( IZylka et al.l 1 19981 ). We reached a la sensitivity of 0.63 mJy and the 
source is detected at 8.29 mJy(^ 13cr), making it the brightest 250 GHz detection among 
all the known z^6 quasars. 



We then searched for the CO (6-5) line emission from this object using the PdBI. The 
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observation was carried out in 2011 in D configuration with a synthesized beam size (FWHM) 
of 5.4^' X 3.9'^ We spent 3 hours on source using the new wide band correlator WideX, with a 
bandwidth of 3.6 GHz in dual polarization. The phase was checked with the quasar 3C 454.3 
and the fiux calibration uncertainty is 5%. The data were reduced with the IRAM GILD AS 
software package (Guilloteau & Lucas 2000). We detect the line at 11a with a typical rms of 
0.5 mJy Beam" ^ per 100 kms~^ channel (Table 1 and Figure 1). We also detect the 99GHz 
dust continuum at 8a with a continuum la sensitivity of 0.05mJy Beam~^. 



2.2. ALMA observations 

We selected J2310+1855 and another four z^6 q uasars that were detected in bright 



millim eter dust continuum and CO (6-5) line emission ( jPetric et al.ll2003l : IWang et al.ll2008L 



2010 : Fan et al. 2013 in prep.) and have Declination of 5 < +20°, so they can be observed 
with ALMA. We list the sources and measurements from previous millimeter observations 
in Table 1. The ALMA observations were carried out in Cycle in 2012 with 14 to 21 12-m 
diameter antennas in the extended array. The typical synthesized beam size (FWHM) is 
corresponding to 4 kpc at z^6. We observed the [C II] 158 /xm (1900.5369 GHz) line 
emission in Band 6 or 7, with one 2 GHz window centered at the line frequency and the 
other three 2 GHz windows observing the dust continuum. The correlator channel width is 
15.625 MHz, or 16^18 kms~^ at the observed line frequency. We spent 50 to 90 minutes 
on source for each target and reached la spectral sensitivities of 0.4 to 0.7 mJyBeam"^ per 
62.5 MHz re-binned channel after continuum subtraction. The continuum sensitivities are 
0.05 to 0.09 mJyBeam-^ 



We reduced the data using the CAS A 3.4.0 packag4j. The fiux scale was calibrated in 
four of the sources by observing Titan, Neptune, and Uranus, and the calibration uncer- 
tainties are better than 15%. For the source J1044— 0125, no primary amplitude calibrator 
was observed. The fiux scale was determined with 3C 273 based on the ALMA and SMA 
observations in March and April 2012 at 225 GHz and 343 GHz. Considering the uncer- 
tainties in the fiux densities and spectral index of the calibrator, we estimate a total fiux 
calibration error of 20% for the observation of J1044— 0125. We derive the line frequencies in 
the kinematic local standard of rest (LSRK) frame to compare the results to previous PdBI 
CO (6-5) line observations. 



http:/ /casa. nrao.edu/ 
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3. Results 

We detected strong [C II] line and dust continuum emission in the host galaxies of all 
the five z^6 quasars. The line emission of the five objects and the continuum emission from 
four of them (the exception is J1044— 0125) are marginally resolved, suggesting deconvolved 
full width at half maximum (FWHM) major axis sizes of 0.3^' to 0.6^' (1.7 to 3.5 kpc) for the 
line and 0.2^' to 0.4^' (1.2 to 2.3 kpc) for the dust continuum. The line intensity-weighted 
velocity maps of J0129— 0035, J1319+0950, and J2310+1855 show clear velocity gradients 
(Figure 2). We subtract the underlying dust continuum in the U— V plane, calculate the line 
spectra by integrating the intensity over the line-emitting region in each channel, and derive 
the redshift, FWHM line width, and line flux by fitting a Gaussian profile to the spectra. 
The [C II] results are summarized in Table 2, and the line and continuum maps are presented 
in Figure 2, and 3. We describe the detailed results for each target below. 

SDSS J2310+1855. This object is the brightest detection among all the z^6 quasars 
from our MAMBO observations. The 250 GHz continuum fiux density is 8.29 ± 0.63 mJy. 
Our PdBI observation of the CO(6-5) line emission shows a line fiux of 1.52±0.13 Jykms"^, 
which is also by far the strongest CO line fiux. The redshift and FWHM measured with 
the CO (6-5) line are zco = 6.0025 ± 0.0007 and FWHMco = 456 ± 64kms-^ The dust 
continuum at 99 GHz has also been detected in the line-free channels with a fiux density of 
0.40±0.05 mJy (Figure 1). 

We detect the [C II] line emission with ALMA in 50 minutes on-source time with a fiux 
which is again more than twice stronger than the four other sources. We fit the spatial extent 
of the line emission in the velocity-integrated intensity map with a 2-dimensional Gaussian 
profile using the IMFIT task in CASA. The derived FWHM source size, after deconvolving 
the 0.73'' X 0.52'' synthesized beam, is (0.56" ± 0.03") x (0.39" ± 0.04") with a position angle 
(PA) of 142° ± 179° east of north. The line intensity- weighted velocity map derived with the 
line-emitting channels indicates a velocity gradient from Southwest to Northeast (See the 
right panels of Figure 2). The 263 GHz continuum fiux density imaged with the line- free 
windows and fitted to a 2-D Gaussian is 9.1 ± 1.4 mJy, which is in very good agreement 
with the MAMBO value. This constrains the (deconvolved) continuum source size to be 
(0.25" ± 0.01") X (0.20" ± 0.02") with PA = 164° ± 178°. 

J1319+0950. This quasar was discovered in the UKIRT Infrared Deep Sky Survey 



(UKIDSS) and has iru^q = 19.65 fl Mortlock et all 120091 ). We detect the [C II] line emission 



in 80 minutes on-source time. The synthesized beam of the line intensity map is 0.69" x 0.49" 
and a 2-D Gaussian fit suggests line emission with a deconvolved source size of (0.58" ± 
0.07") X (0.32" ± 0.16") oriented at PA = 27° ± 18°. The line velocity map shows a clear 
gradient from the Southwest to the Northeast. We plot the [C II] line channel map of this 
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object in Figure 3. The map shows a clear position shift of the [C II] hne peak (by ~ 0.4'^) 
from the +222kms~^ channel to the — 200kms~^ channel. The 258 GHz continuum flux 
density is 5.2 ± 0.8 mJy and the deconvolved continuum source size is (0.37^^ ± 0.02'^) x 
(0.32^^ ± O.O3'0 with PA = 102° ± 176°. 

J2054— 0005. The quasar was selected from SDSS stripe 82 with mi45o = 20.60, 
i .e. about one ma gnitude fainter than the objects discovered from the SDSS main survey 
( j Jiang et al.l 120081 ). We observed the [C II] line emission for 40 minutes. The line emission 
is marginally resolved with the 0.65^' x 0.59^' synthesized beam, with a deconvolved FWHM 
source size of (0.35'' ± 0.05'') x (0.32" ± 0.06") and PA = 82° ± 54° The 262 GHz contin- 
uum flux density measured in the line- free windows is 3.1 ± 0.05 mJy, and the deconvolved 
continuum source size is (0.28" ± 0.03") x (0.26" ± 0.03") with PA = 168° ± 171°. 

J0129— 0035 . This quasar was discovered in SDSS Stripe 82 with mi45o = 22.16 



( 1 Jiang et al.l 120091 ). It is the faintest optical source among the MAMBO-detected z^6 
quasars, but previous millimeter obse rvations found brig ht 250 GHz dust continuum and CO 
(6-5) line emission from this object ( jWang et al.ll2011al ). The [C II] line was detected with 
1-hour on-source time with ALMA. The source size fltted with a 2-D Gaussian proflle and 
deconvolved by a 0.57" x 0.49" beam is (0.41" ±0.06") x (0.25" ±0.10") with PA = 11° ±161°. 
A velocity gradient is found from the Southeast to the Northwest. A 2-D Gaussian flt to 
the continuum image gives a 287 GHz flux density of 2.58 ± 0.05 mJy, with a deconvolved 
continuum source size of (0.23" ± 0.03") x (0.14" ± 0.06") at PA = 115° ± 179°. 

J1044— 0125. This is another broa d absorption lin e quasar at z^6 discovered from the 
SDSS main survey with mi45o = 19.2 (IFan et al.l l2000l ). A deep near-IR spectrum of the 
quasar C IV line emission of this object yields a SMBH mass of 10.5 x 10^ Mq (Jiang et 
al. 2007). We detected the [C II] line in 87 minutes of on-source time. The beam size of 
the [C II] line observation is 0.66" x 0.45" and the deconvolved source size is estimated to 
be (0.61" ± 0.11") X (0.33" ± 0.22") oriented at PA = 65° ± 48°. No gradient is seen in the 
line intensity- weighted velocity map of this object. The 287 GHz continuum flux density is 
3.1 ± 0.6 mJy and the source is unresolved in the continuum. 

We c alculate the [C II] line luminosity of each source asL[cn]/L0 = 1.04x 10~^SAvi/o(l + 
z)-iDl^ J Solomon fc Vanden Boutlboos h. where z/q = 1900.5369 GHz is the rest frame [C 
II] line frequency, SAv is the integrated line flux in Jykms"^, and Dl is the luminosity 
distance in Mpc. The derived line luminosities are in the range of 1.6 to 8.8xlO^L0. We 
also calculate the CO (6-5) line luminosity for J2310+1855, which is (5.4 ± 0.5) x 10^ Lq or 
(5.1 ± 0.4) X 10^^ Kkms~^ pc^. We then estimate the FIR luminosities by fltting a modifled 
black body to the ALMA continuum flux densities, the MAMBO 250 GHz flux densities, 
and available PdBI 100 GHz and SCUBA 350 GHz data. As there are no measurements at 
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shorter wavelengths to constrain the peak of the FIR SED and thereby the dust temperature, 
we adopt a dust temperature of 47 K and emissivity in dex of 1.6, which ar e the mean values 
found in the sample of high-z FIR luminous quasars (IBeelen et al.l l2006l ) . This yields the 
FIR luminosities in the range 42.5-122.5 /xm of 5— 17x10^^ L©. We have listed the derived 
[C II] and FIR luminosities in Table 3. 



4. Discussion 

4.1. Distributions of gas, dust, and star formation in the quasar host galaxies 



The new detections of strong [C II] line emission toward the five FIR-luminous quasars 
at z^6 provide further evidence for active star formation in the quasar host galaxies. The [C 
II] line emission is marginally resolved in all these objects. The ALMA observations at 0.7^' 
resolution yield estimates of the intrinsic [C II] source sizes (FWHM of the major axis) of 
0.3^' to 0.6^' (1.7 to 3.5 kpc). Clear velocity gradients have been found in the line intensity- 
weighted velocity maps of J2310+1855, J1319+0950, and J0129— 0035, which suggests that 
the gas could be in rotation, and gravitationally bound, e.g. a non-face-on rotating disk. The 
dust continuum fiux densities measured in the vicinity of the [C II] line frequencies for all five 
objects are consistent with previous MAMBO measurements at 250 GHz. The continuum 
sources of four of them (except J1044— 0125) are marginally resolved, indicating deconvolved 
FWHM major axis sizes of 0.2^'' to 0.4'^, or 1.2 to 2.3 kpc. These results constrain the spatial 
extent of star forming activity to be 2.6 to 5.3 kpc in diameteiil in the nuclear region. We 
will observe these sources with ALMA in Cycle 1 at 0.2'^ resolution to fully resolve the line 
and dust continuum sources, and measure what fraction of the dust continuum emission is 
from the central compact AGN. These observations will finally measure the star formation 
rates (SFR) and SFR surface densities in these earliest quasar-starburst systems. 

In Figure 4, we plot the [C II]-to-FIR luminosity ratios of of the [C Ilj-detected z>5. 8 

quasars, including the five new detections in this work, J1148+5251 (IMaiolino et al.ll2005h. 

CFH QS J0210-0456 at z=6.43 flwillott et al.l[2013l ). and ULAS J1120+0640 at z=7.08 dVenemans e^ 
20121). We compare them to samples of [C II]-detected local normal star f orming galaxies, 



ULIRGs, submillimeter galaxies, and FIR-luminous quasars at high redshift (IMalhotra et al 



2001 



Luhman et al.ll2003: Stacev et al 



2010: Maiolino et al.ll2009 



Ivison et al 



2010 



De Breuck et al. 



2011 



Swinbank et al 



2012 



Wagg et al. 



20121 : 



Petv et al.l[20M : iGallerani et aPboid : ICariUi et al. 



^We adopt a source size of 1.5 x the FWHM major axis from the [C II] intensity map (i.e., full width at 
20% of the peak intensity for a Gaussian profile). 
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20131 ). The five z^6 quasars presented in tliis work, as well as J 1148+5251, show luminosity 
ratios of 3 to 5.6 x 10""^, which are comparable to the typical values found in local ULIRGs 
and 1 < z < 5 [C II]-detected quasars, and a few to 10 times lower than that of the disk star 
forming galaxies and submillimeter galaxies. We also notice that the other two z>6 quasars 
with moderate FI R luminosities (10-^ -^ to 10^^ L©) show higher L[cn]/LFiR than most FIR 
luminous objects (jWillott et al.ll2013l ). We estimate the CO (1-0) line luminosities for the 
five objects from the PdBI CO (6-5) det ections (last Column of Table 1), assuming a CO 
excitation ladder similar to J1148+5251 (IRiechers et al.ll2009l ). The calculated [C II]-to-CO 
(1-0) line luminosity rat ios are about 2500 to 4700, which is slightly higher than the values 
found in local ULIRGs (ILuhman et al.ll2003l ). and the highest value is close to the median 



luminosity ratios found in starburst galaxies (i.e. 
2010l : ISwinbank et al.l[2012l ). 



L[cn]/Lco ^ 4400, IStacey et al 



1991 



The [C II], FIR, and CO luminosity ratios of thes e [C II]-detected z^6 



be produced by photodiss ociation region (PDR) models (IKaufman et al.lll999l : IStacey et al. 



quasars can 



2010l : iLuhman et al.ll2003l ) with gas density on orders of 10^ to 10^ cm ^ and FUV (6 eV< 
hu < 13.6 eV) rad iation field Go of a fe w 10^ to 10^ (in units of the Habing Field, 1.6 x 



10~^ergscm~^s~^, iKaufman et al.lll999l ). However, such a strong radiation field will also 
prod uce strong [O I] 63/im line emission with intensity comparab le or higher than the [C II] 
line (IKaufman et al.lll999l : lLuhman et al.ll2003l : ICarilli et al.ll2013l ). Thus future observations 



of the [O I] and other fine structure lines from these objects with the full configuration of 
ALMA will be a crucial test of the physical conditions inferred from these PDR models. It 
is also possible that the central AGN has substantial contribution to the dust heating and 
FIR emission, which results in lower [C I I]-to-FIR luminosity ratios in the nuclear region 
(ILuhman et al.ll2003l : ISargsyan et al.ll2012l ). 



We compare the [C II] and CO (6-5) line profiles of the five objects in Figure 5. For four 
of the five sources, the redshifts measured with [C II] and CO (6-5) are consistent within the 
la errors; there are no large velocity offsets between the [C II]-traced atomic gas and the CO 
(6-5)-traced molecular gas. The [C II] FWHM line widths of J2310+1855, J2054-0005, and 
J0129— 0035 are about 60 to 115 kms~^ smaller than the CO (6-5) measurements. But these 
differences are within the 1 to 2 a error bars, as the CO (6-5) line width uncertainties for 
these objects are between 60 and 120 kms~^ (Table 1). The other object J1044— 0125, shows 
a larger [C II] redshift with Az = Z[cii] — zco = 0.0023 ± 0.0010 (i.e. a velocity difference 
of 100 ± 44kms~^), and a much broader [C II] line width (Figure 5). This may indicate 
different kinetic properties between the atomic and molecular gas components in this object. 
However, it is also possible that a large fraction of the CO line emission is undetected and 
the CO line width is underestimated due to the low signal-to-noise ratio of the line spectrum. 
Thus, deep imaging of the CO line emission with better measurements of the CO line profile 
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and spatial distribution is needed to address whether the [C II]-emitting gas is more centraUy 
concentrated in J2310+1855, J2054— 0005, and J0129— 0035, and to understand the origin 
of the narrower CO hne in J1044— 0125. 



4.2. Dynamical masses traced by the [C II] line emission 



If we assume a rotating disk geometry for the [C II] -emitting gas in these FIR-luminous 
z^6 quasar^, we can estimate the dynamical masses within the [C II] -emitting regions for 
J2310+1855, J1319+0955, J2054-0005, and J0129-0035 as Mdyn/M© ^ 1.16 x lO^Vdr^D, 
where D is the disk diameter in kpc from the [C II] measurements (2.6 to 5.3 kpc) and 
Vcir is the maximum circular velocity of the gas disk in km s~^. We estimate Vcir as Vcir = 
0.75FWHM[cii]/sinz (i.e., half width at 20% line maximum), where i is the inclination angle 
between the gas disk and the line of sight (i = 0° for a face-on disk). The derived M^ynSin^ i 
are 0.9 — 8.8 x 10^° M©, and the uncertainties estimated with the measurement errors in 
FWHM[cn] and [C II] source size are about 10% to 40%. With the assumption of an inclined 
disk geometry, we can also have a first guess at the disk inclination angle from the [C 
II] minor and major axis ratios (R), i.e. i = cos~-'^(R), though we should keep in mind 
that the source size measurements at current spatial resolution (^ 0.7^') still have large 
uncertainties (see Section 3). The estimated inclination angles are 46°, 57°, 24°, and 52° for 
J2310+1855, J1319+0955, J2054-0005, and J0129-0035, respectively. This gives inclination 
angle-corrected dynamical masses of M^yn = 1.4 x 10^° to 1.3 x 10^^ M©. We list the derived 
inclination angle-corrected dynamical mass M^yn in Table 3. 

There are no published SMBH masses for the four objects yet. We calculate the total 
quasar bolometric luminosities for the four sources using an isotropic bolometric correction 
of Lboi = 4.2Li45o ( IRunnoe et al.ll2012l ). where L1450 is the rest-frame 1450 A luminosity, and 
derive the SMBH masses (Mbh) from the quasar luminosity assuming Eddington accretion. 
The resulting mass ratios Mbr/M^j^ lie in the range of 0.012 to 0.030 (see Table 3). Note that 
the Mbh and Men/Mdyn values could be even larger if the typical rat io between quasar bolo - 
metric luminosity and Eddington luminosity is less than unity (e.g. iDe Breuck et al.ll201ll ). 
This is consistent with the previous CO estimates of the median SMBH-bulge mass ratio of 
these z^6 quasars and agrees with the results found with other high-z FIR 

and CO luminous quasars that the SMBH-bulge mass rat ios are 10 to 30 times higher than 



the average value of 0.0014 found in local normal galaxies (jMarconi et al.ll2003l : IWalter et al 



^We exclude J1044— 0125 in the discussion here as the difference between the [C II] and CO (6-5) line 
spectra may indicate more comphcated gas dynamics. 
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20041 : 



Riechers et al 



2008 



Wang et al.ll2010l : ICoppin et al. 



2008 



Venemans et al.ll2012l ). 



However, one should be cautious with these [C II]-based Mbr/M^j^ values, as the de- 
tected [C II] line emission may trace only the intense star forming region in the very centers 
of quasar host galaxies and not extend as far as the stellar bulge. This may underestimate 
the total masses within the spheroidal stellar bulges. Additionally, we notice that there are 
still large uncertainties in the intrinsic source morphology measurements; the deconvolved 
minor axis and position angle measurements of some of the objects show large error bars. 
This can introduce significant uncertainties in the inclination and dynamical mass estimate. 
Further ALMA imaging of these FIR-luminous z^6 quasars at ^ 0.2^' will be crucial to better 
constrain the distribution and dynamical properties of the gas components, address how well 
the atomic/molecular line emission traces the dynamical masses of the quasar host galaxies, 
and determine whether the high MsH/Mdyn ratios is common in the massive quasar-starburst 
systems in the early universe. 



4.3. Search for [CII] line emitters associated with the quasar environments 

Given the sensitivity of the ALMA Band 6 and Band 7 data, and the large spectral 
bandwidth (7.5 GHz, or ^9000 kms~^), it is plausible that serendipitous line emission from 
companion objects will be detected in our data. Although the volume is relatively small 
over the five fields (each field was observed across four 1.875GHz spectral windows, cover- 
ing a redshift range Az ^ 0.2 which, for a primary beam with FWHM^ 22^' results in a 
total surveyed volume of ^ 280 Mpc^), recent ALMA surveys of submm luminous starburst 
galaxies at z ^ 4.4 suggest a strong evolution in the [C II] luminosity function out to these 
early cosmic times (Swinbank et al. 2012). We use a Bayesian search algorithm developed 
to eflSciently search for broad, weak line emission in large spectral line data cubes (Lentati 
et al. 2012), and use the evidence for candidate detections to calculate the probability that 
detected candidates are 'real' relative to being noise. We do not find any candidate detec- 
tions of serendipitous line emission in our data with probabilities above 25%. At the current 
sensitivity the volume sampled is not sufficient to place constraints on the possible evolution 
of the luminosity function out to z ^ 6, with source counts (i.e., 10~^Mpc~^, Swinbank et 
al. 2012) derived from the [CII] luminosity function predicting ^ 0.3 detection above a 3a 
luminosity limit L[cn] > 2.1 x 10^ Lq. However the sensitivity of ALMA at band 6 and 7 
means that deeper observations in Cycle 1 of similar volumes will allow us to place strong 
constraints on the obscured star-formation properties of galaxies in the environments of the 
quasars. 
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5. Summary 

We detected [C II] fine structure line and dust continuum emission from the host galaxies 
of five quasars at z~6, using ALMA in the Cycle phase at ~ 0.7^' resolution. Complemen- 
tary IRAM observations of the CO (6-5) line and 250 GHz dust continuum emission from 
the z=6.00 quasar J2310+1855 are also presented. Our ALMA observations measure the 
FWHM major axis sizes of the [C II] emission from the five objects to be 0.3'' to 0.6'' (1.7 to 
3.5 kpc) and the sizes of the dust continuum source for four of them to 0.2" to 0.4" (1.2 to 
2.3 kpc). The detections of [C II] line and dust continuum emission suggestiong star forming 
activity in the central few kpc region of the quasar host galaxies. The derived [C II]-to-FIR 
luminosity ratios are of the order 10~^, which are comparable to the typical values found 
in local ULIRGs and other FIR-luminous quasars at high redshift. The intensity- weighted 
velocity maps of J2310+1855, J1319+0950, and J0129-0035 show velocity gradients. Such 
velocity gradients are consistent with rotation, although they are not uniquely interpreted as 
such. We estimate the dynamical masses within the [C II]-emitting region for four of the five 
objects assuming that the gas is distributed in a rotating disk. The derived ratios between 
the SMBH masses and the dynamical masses are one order of magnitude higher than that 
of local normal galaxies. 

The detections of [C II] 158/im line emission from quasar host galaxies at z^6 have 
demonstrated the power of ALMA in observing signatures of star formation at the earliest 
cosmic epoch. With the full configuration of ALMA, we should be able to resolve the hne and 
dust continuum emission on kpc or sub-kpc scales, which will measure the surface densities 
of the atomic/molecular gas components and star forming activity in the nuclear region, and 
better address the gas kinetics and dynamical masses of the spheroidal quasar host galaxies. 
The full frequency coverage of ALMA will also allow a search of other ionized/atomic inter- 
stellar coohng lines from these FIR-luminous z^6 quasars to measure the physical conditions 
(e.g. density, rediation field, temperature, etc.) of the interstellar medium in these earliest 
quasar-starburst systems. 

We thank M. Lacy at the National Radio Astronomy Observatory for help with the 
observation and data analysis. This work is based on observations carried out with ALMA 
(NRAO), the Max Planck MiUimeter Bolometer Array (MAMBO) on the IRAM 30m tele- 
scope, and the Plateau de Bure Interferometer. The National Radio Astronomy Observatory 
(NRAO) is a facility of the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. This paper makes use of the following ALMA data: 
ADS/JAO.ALMA# 2011.0.00206. S . ALMA is a partnership of ESO (representing its mem- 
ber states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and 
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ASIAA (Taiwan), in cooperation with the Repubhc of Chile. The Joint ALMA Observatory 
is operated by ESO, AUI/NRAO and NAOJ. IRAM is supported by INSU/CNRS (France), 
MPG (Germany) and IGN (Spain). Frank Bertoldi and Fabian Walter acknowledge support 
through the DFG priority program 1573 and the SFB 956. Desika Narayanan acknowledges 
support from the NSF via grant AST- 1009452. Kirsten Knudsen acknowledges support from 
the Swedish Research Council. X. Fan acknowledges support from NSF grant AST 08-06861 
and 11-07682 and a David and Lucile Packard Fellowship. 
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23^0"'40' 39' 38' 23^ 0"'40' 39' 38' 



Fig. 1. — PdBI observation of CO (6-5) line and dust continuum emission at 99GHz from 
J2310+1855. The left panel shows the CO (6-5) line spectrum binned to 30 kms~^ channels. 
The solid line is a Gaussian fit to the line spectrum. The middle panel shows the intensity 
map of the CO (6-5) line emission. The la rms noise of the map is 0.13 Jykms"-^ Beam~^ 
and the contours in steps of 2a. The beam size of 5.4^' x 3.9^' is plotted on the bottom left. 
The cross denotes the position of the optical quasar. The right panel is the continuum map 
at 99 GHz with la rms of 0.05mJy Beam~^, and the contours in steps of 2a. 
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Table 1: Summary of previous observations 



Source 


mi450 


S250 


zco 


FWHMco 


SAvco(6-5) 


Lco(i-o) 






mJy 




kms"-*^ 


Jy kms"-*^ 


105 Lo 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSS J231038.88+185519.7 


19.30 


8.29±0.63 


6.0025+0.0007 


456+64 


1.52 + 0.13 


32.1 


ULAS J131911. 29+095051. 4 


19.65 


4.20±0.65 


6.1321+0.0012 


537+123 


0.43 + 0.09 


9.4 


SDSS J205406.49-000514.8 


20.60 


2.38±0.53 


6.0379+0.0022 


360+110 


0.34 + 0.07 


7.3 


SDSS J012958.51-003539.7 


22.28 


2.37±0.49 


5.7794+0.0008 


283+87 


0.37 + 0.07 


7.4 


SDSS J104433. 04-012502. 2 


19.21 


1.82±0.43 


5.7824+0.0007 


160+60 


0.21 + 0.04 


4.2 



Note - Column (1), name; Column (2), magnitudes at rest-frame 1450 A (Fan et al. 2013. in prep: iMortlock et al.ll2009l : 
I Jiang et al.ll2008l . I2OO9I : IPan et al.ll2000 h: Column (3), continuum flux density at 250 GHz; Column (4) and (5), redshift and 
FWHM line width of the CO (6-5) hue; Column (6) CO (6-5) hue flux of the flve objects; Column (7) CO (1-0) line luminosity 
calculated from the CO (6-5) line emission assuming a CO excitation ladder similar to J1148+5251 jRiechers et al I I2OO9I ). The 
CO (6-5) line measurements for J2310+1855 are from this work, and from lWang et al.l (l2010l . l2011al ) for the other four sources. 



Table 2: ALMA observations 



Source 


ton 


Z[CII] 


SAv 


FWHM[cii] 


rms[cii] 


^con 


Scon 


rmscon 




min 


Jy kms"-*^ 


km s"-*^ 


mJy Beam"-*^ 


GHz 


mJy 


mJy Beam"-*^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


J2310+1855 


50 


6.0032+0.0002 


8.95+1.42 


392+22 


0.5 


263 


9.1+1.4 


0.06 


J1319+0950 


80 


6.1329+0.0007 


4.27+0.88 


518+81 


0.7 


258 


5.2+0.8 


0.08 


J2054-0005 


40 


6.0391+0.0001 


3.45+0.53 


242+10 


0.4 


262 


3.1+0.5 


0.04 


J0129-0035 


60 


5.7787+0.0001 


2.00+0.32 


194+12 


0.4 


287 


2.6+0.4 


0.05 


J1044-0125 


87 


5.7847+0.0007 


1.70+0.45 


420+75 


0.6 


287 


3.1+0.6 


0.09 



Note - Column (1), source; Column (2), on-source time of the ALMA [C H] observations; Column (3), (4), and (5), redshift, 
flux, and FWHM line width of the [C H] line fltted to a single Gaussian line proflle; Column (6), ALMA la line sensitivity bin 
to a channel width of 62.5 MHz (~ 70kms~-'^); Column (7) and (8), frequency and flux density of the continuum measured 
with the line- free windows; Column (9), ALMA la rms continuum sensitivity integrated over a total bandwidth of 5.8 — 6 GHz. 
The 15% to 20% calibration uncertainties are included in the flux errors. 
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Table 3: Luminosities and Dynamical masses 



Source 


L[cii] 


Lfir 


Lbol 


Mbh 


MdynSin^i 


Mdyn 


MBH/Mdyn 




109 L0 


1012 L0 


1013 L0 




IOIOM0 


1O1OM0 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


J2310+1855 


8.8±1.4 


17.2±1.0 


9.3 


2.8 


4.9±0.6 


9.6 


0.030 


J1319+0950 


4.3±0.9 


10.4±1.0 


7.0 


2.1 


8.8±3.0 


12.6 


0.017 


J2054-0005 


3.4±0.5 


6.1±0.8 


2.8 


0.86 


1.2±0.2 


7.1 


0.012 


J0129-0035 


1.9±0.3 


4.7±0.6 


0.57 


0.17 


0.9±0.2 


1.4 


0.012 


J1044-0125 


1.6±0.4 


5.3±0.6 


11.6 


10.5 









Note - Column (1), source; Column (2), [C II] line luminosity; Column (3), FIR luminosity in the wavelength range from 
42.5/im to 122.5/im. Column (4), quasar bolometric luminosity calculated from the rest-frame 1450 A magnitude, assuming 
an isotropic bolometric correction of Lbol = 4.2Li45o (jRunnoe et al.l [2^0121 ): Column (5), SMBH masses. The SMBH mass 
of J1044— 0125 is calculated from the quasar C IV line emission (Jiang et al. 2007), and we estimate the SMBH masses for 
the other four sources from the bolometric luminosities assuming Eddington accretion; Column (6), dynamical mass without 
inclination angle correction, estimated from the [C II] line width and source size (see Section 4.2); Column (7), inclination 
angle-corrected dynamical mass within the [C II] -emitting region. We estimate the disk inclination angle from the [C II] minor 
and major axis ratio as i = 46°, 57°, 24°, and 52° for J2310+1855, J1319+0950, J2054-0005, and J0129-0035, respectively; 
Column (8), SMBH-bulge mass ratio. 
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13"l9"'ir.45 ir.35 ir.25 ir.15 
J2000 Right Ascension 




20^54"'06^65 06^55 06^45 06^35 
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Fie. Z. The dust continuum map (left), [C II] line velocity-integrated map (middle), and line velocity maps (right) of the 

five new [C II] detections. We calculate the line intensity- weighted velocity map using pixels detected at > 4cr in each case. 
The black crosses show the position of the optical quasar from the discovery paper (Fan et al. 2013, in prep., [Mort lock et al.l 
[2009; Jiang et al. 2008, 2009; Fan et al. 2000). The sizes of the synthesized beams are plotted in the bottom-left of each panel. 
J2310+1855 - The continuum contours are [-2, 2, 4, 8, 32, 64] x 0.1 mJy Beam" \ and the line contours are [-2, 2, 4, 8, 16, 
32]x0.15 Jy Beam~^ kms~^. The velocity contours are [-1, 0, 1, 2, 3]x40 km where zero velocity corresponds to the CO 
(6-5) redshift of zco = 6.0025. The la rms noise is 0.06 mJy Beam" -"^ for the continuum map and 0.14 Jy km s"-*^ Beam" -"^ 
for the line map. J1319+0950 - The contours are [-2, 2, 4, 8, 16, 32] x 0.1 mJy Beam" ^ for the continuum, [-2, 2, 4, 6, 8, 
10, 12] X 0.18 J y Beam~^ kms~^ for the line, and [-2, -1, 0, 1, 2, 3]x75 km for the velocity map. Zero velocity is defined 
at Zco = 6.1321. The la rms noise values are 0.08 mjy Beam" -"^ and 0.18 Jy km s"-*^ Beam" -"^ for the continuum and line 
maps, respectively. J2054-0005 - The line contours are [-2, 2, 2.83, 4, 5.66, 8, 11.31, 16, 22.63] x0.15 Jy Beam"^ /cm5"\ 
the continuum contours are [-2, 2, 4, 8, 16, 32, GAjxOMmJy Beam~^ , and the contours in the velocity map are [0, 1, 
2]x4:0 km . Zero velocity corresponds to zco = 6.0379. The la rms noise is 0.04 mJy Beam" -"^ for the continuum map 
and 0.08 Jykms"^ Beam"^ for the line map. J0129-0035 - The contours are [-2, 2, 4, 8, 16, S2]x0.1 mJy Beam~^ for the 
continuum, [-2, 2, 4, 8, 16, 32] x 0.075 Jy Beam"^ /cms"^ for the line, and [-3, -2, -1, 0]x20/cms"i for the velocity map. The 
zero velocity is set to zco = 5.7794. The la rms noise is 0.05 mjy Beam" -"^ for the continuum map and 0.08 Jykms"^ Beam"^ 
for the line map. J1044— 0125 - The contours are [-2, 2, 4, 8, 16] x 0.08 mJy Beam" -"^ for the continuum and [-2, 2, 4, 6, 
8] X 0.14 Jy Beam" ■'^ kms~^ for the line. Zero velocity is set to zco = 5.7824. The la rms noise is 0.09 mjy Beam" -"^ for the 
continuum map and 0.14 Jy km s"-*^ Beam"-*^ for the line map. 



- 21 - 




Fig. 2. — Continued 
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Fig. 3.^ [C II] line channel maps of J1319+0950 in the velocity range of -340.3 to +362.9 km s"^ at 
0.7^' X 0.5^' resolution. The channel width is 140kms~^ for each map, and the line is clearly detected 
in the middle four channels. The typical la rms noise of the channel map is 0.5 mJy Beam~^, and 
the contours are [-2, 2, 4, 6, 8, 10, 12]xlcr. The peak of the [C II] line emission is shifted by about 
0.4^' from the 222kms~^ channel to the — 200kms~^ channel. The cross denotes the location of 
the optical quasar. 
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Fig. 4. — L[cii]/Lfir vs. Lfir. The black stars indicate the quasa rs at z>5.8, includin g the five 
new [C II]-detected z~6 quasars in this work , J1148+5251 at 7^6.42 (lMaio1ino et a,l.ll2nn5l v CFHQS 
J0210-0456 at z=6.43 (Iwillott et aklboLsh . and the z=7.08 quasar ULAS J112001.48+064124.3 
(jVenemans et al.l 120121'). We also plot the lumino sity ratios from sarn ples of local normal star 
forming gala xies dMalhotra et al.l [2OO1I ). ULIRGs f Luhman et al. 2003 1, 1 < z < 5 star form- 
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Fig. 5. — The [C II] line spectra (black solid line) of the five ALMA detected quasars integrated 
over the line-emitting region in each channel, together with the previous CO (6-5) detections from 
PdBI (blue dotted hne, scaled to the [C II] hne) and a Gaussian fit to the [C II] line (red dashed 
line). The typical ±lcr error per channel for the [C II] line spectrum is shown in the top right of 
each panel. 



